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We discuss a possible origin of the experimentally observed nonlinear contribution to the shift 
ATc = Tc — T® of the critical temperature Tc in an atomic Bose-Einstein condensate (BEC) with 
respect to the critical temperature T® of an ideal gas. We found that accounting for a nonlinear 
(quadratic) Zeeman effect (with applied magnetic field closely matching a Feshbach resonance held 
Bo) in the mean-held approximation results in a rather signihcant renormalization of the held-free 
nonlinear contribution 62 , namely ATc/T° ~ fe|(a/AT)^ (where a is the s-wave scattering length, At 
is the thermal wavelength at T°) with 62 — 7^52 and 7 = 7 (Bo). In particular, we predict 62 — 42.3 
for the Bo — 403G resonance observed in the BEC. 

PACS numbers: 67.85.Hj, 67.85.Jk 


Studies of Bose-Einstein condensates (BECs) continue 
to be an iimortant subject in modern physics (see, 
e.g., Refs.[ll-y] and further references therein). Atomic 
BECs are produced in the laboratory in laser-cooled, 
magnetically-trapped ultra-cold bosonic clouds of differ¬ 
ent atomic species (including a ‘^^Na 

0, 0 4^6 O, ” 

and [17 


41 


K [ill, 433(7-3 [121, [13| 

among others). Also, a discussion of a 
relativistic BEC has appeared in Ref.[l^ and BECs of 
photons are most recently under investigation [l^ . In ad¬ 
dition, BECs are successfully utilized in cosmology and 
astrophysics [I3 as they have been shown to constrain 
quantum gravity models [i0 . 

In the context of atomic BECs interparticle interac¬ 
tions must play a fundamental role since they are nec¬ 
essary to drive the atomic cloud to thermal equilibrium. 
Thus, they must be carefully taken into account when 
studying the properties of the condensate. Eor instance, 
interatomic interactions change the condensation tem¬ 
perature Tr, of a BEC, as was pointed out first by Lee 
and Yang [ii [13 (see also Refs.HIHa^ for more recen'i 
works). 

The first studies of interactions effects were focused 
on uniform BECs. Here, interactions are irrelevant in 
the mean field (ME) approximation (see Refs. 25, 
but they produce a shift in the condensation temperature 
of uniform BECs with respect to the ideal noninteract¬ 
ing case, which is due to quantum correlations between 
bosons near the critical point. This effect has been fi¬ 
nally quantified in 25, as ATc/T° ~ 1.8(a/AT), where 
ATc = Tc — T^ with T). the critical temperature of the gas 
of interacting bosons, T° is the BEC condensation tem¬ 
perature in the ideal noninteracting case, a is the s-wave 
scattering length used to represent interparticle interac¬ 


tions 000, and At = is the thermal 

wavelength for temperature with ma the atomic mass. 

But laboratory condensates are not uniform BECs 
since they are produced in atomic clouds confined in mag¬ 
netic traps. For trapped BECs, interactions affect the 
condensation temperature even in the MF approxima¬ 
tion, and the shift in Tc in terms of the s-wave scattering 
length a is given by 


ATc/T^ -bi{a/XT) + b2ia/XTf 


( 1 ) 


with bi ~ —3.4 0 and 62 — 18.8 isil . 

High precision measurements |32| of the condensation 
temperature of 3® AT in the range of parameters A^ ~ (2 — 
8)x103, 10-3 < o/At < 6x10-2 andTc ~ (180-330)nA: 
have detected second-order (nonlinear) effects in ATc/Tc 
fitted by the expression ATc/Tc = b^^{a/XTY with 
b^^ ~ 46 ± 5. This result has been achieved ex¬ 
ploiting the high-held 403G Feshbach resonance in the 
|F, Tot >= |1,1 > hyperhne (HF) state of a conden¬ 
sate where F = S-\-1 is the total spin of the atom with S 
and / being electron and nuclear spin, respectively, and 
mp is the projection quantum number. Thus, the the¬ 
oretically predicted [3l| quadratic-amplitude coefficient 
62 turned out to be in a rather strong disagreement with 
the available experimental data. There have also been 
some efforts to theoretically estimate the correct value of 
62 in the MF approximation by conside ring anharmonic 
and even temperature-dependent traps 33|, which how¬ 
ever have not been too successful. Therefore one could 
expect that a more realistic prediction of the experimen¬ 
tal value of b^^ should take into account some other so 
far unaccounted effects. 

The main goal of this paper is to show that, taking into 
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account the nonlinear (quadratic) Zeeman effect and us¬ 
ing the MF approximation, it is quite possible to explain 
the experimentally observed [Hj value of 62 for the dOSG 
resonance of the hyperfine \F,mF >= | 1,1 > state of 
with no need to go beyond-MF approximation. 

Recall that experimentally the s-wave scattering length 
parameter a is tuned via the Feshbach-resonance tech¬ 
nique based on Zeeman splitting of bosonic atom levels 
in an applied magnetic field. This means that the in¬ 
teraction constant g = a/ma) is actually always 

field-dependent. More explicitly, according to the inter¬ 
pretation of the Feshbach resonance 3J, 


where abg is a so-called background value of a, Bq is the 
resonance peak field, and A the width of the resonance. 

Thus, in order to properly address the problem of 
condensation-temperature shifts (which are always ob¬ 
served under application of a nonzero magnetic field B), 
one must account for a Zeeman-like contribution. It 
should be emphasized, however, that a single (free) atom 
Zeeman effect (induced by either electronic or nuclear 
spin) g,aB is not important for the problem at hand sim¬ 
ply because it can be accounted for by an appropriate 
modification of the chemical potential. 

Recall that in the presence of a linear Zeeman effect, 
the basic properties of an atomic BEC can be understood 
within the so-called ” condensate wave function” approx¬ 
imation 0 


a clear-cut experimental evidence [H, in favor of the 
so-called Breit-Rabi nonlinear (quadratic) HF-mediated 
Zeeman effect 0 in BEC. We are going to demonstrate 
now how this nonlinear phenomenon (which is not a triv¬ 
ial generalization of the linear Zeeman effect) affects the 
BEC properties (including a feasible condensation tem¬ 
perature shift). Recall that in strong magnetic fields, the 
magnetic-field energy shift of the sublevel tof of an alkali- 
metal-atom ground state can be approximated (with a 
rather good accuracy) by the following expression [s^ 


Enlz = A-hf 


El 

hSvhf 


( 5 ) 


where Ahf = 


1 - 


4ml 


( 21 + 1)2 ) and 6vhf is the so-called 
hyperfine splitting frequency between two ground states. 

Now, by repeating the above-mentioned Bogoliubov’s 
procedure, we obtain a rather nontrivial result for BEC 
modification. Namely, it can be easily verified that HF- 
mediated nonlinear Zeeman effect gives rise to the fol¬ 
lowing two equivalent options for the energy minimiza¬ 
tion (based on the previously defined ground state with 
Ez = gno): (a) Enlz oc g'^nl or (b) Enlz oc Ezgm = 
{g,BB)gnQ. As a matter of fact, the choice between these 
two options is quite simple. We have to choose (b) sim¬ 
ply because (a) introduces the second order interaction 
effects (oc g^) which are neglected in the initial Hamilto¬ 
nian (3). As a result, the high-field nonlinear Zeeman ef¬ 
fect produces the following modification of the local BEC 
energy: 


n = 


/ 


(fxH{x) 


( 3 ) 


SHnlz{x) 


2gnodn{x) + AnFgno 


/ fj^sB 
\h6Fhf 


Sn(x) ( 6 ) 


where H(x) = gn^ — Ezn with n{x) = 'I'''"(x)4'(x) being 
the local density of the condensate (4'(x) is the properly 
defined wave function of macroscopic condensate), and 
Ez = t^sB the Zeeman energy (with /is being the Bohr 
magneton). 

Eollowing Bogoliubov’s recipe [s^, let us consider 
small deviation of the condensate fraction from the 
ground state no (the number) by assuming that n{x) ~ 
nQ-\-6n{x) with 6n{x) <C no. Treating, as usually, no and 
6n(x) independently, we obtain from Eq.(3) that in the 
presence of the linear Zeeman effect gno = Ez (meaning 
that Ez is playing a role of the chemical potential fil) 
and, as a result, the BEC favors the following energy 
minimum: 


Therefore, accounting for nonlinear Zeeman contribu¬ 
tion will directly result in a renormalization of the high- 
held scattering length 


5t: 

a = a 


1 + 2"^hf 


g-pB \ 
hSvhf J 


( 7 ) 


Now, by inverting ([2]) and expanding the resulting B{a) 
dependence into the Taylor series (under the experimen¬ 
tally satished conditions abg a and A <C Ho) 


B{a) ~ Bo 




one obtains 


( 8 ) 


6Ho{x) = 2gnoSn{x) — EzSn{x) = gnoSn{x) (4) 


a* ~ 7 a-I-0(a;,g/a, A/Ho) (9) 


Thus, we come to the conclusion that at low mag¬ 
netic helds (where the linear Zeeman effect is valid ), in 
accordance with the available experimental results |37| . 
there is no any tangible change of the BEC properties 
(including g modihcation). On the other hand, there is 


for an explicit form of the renormalized scattering length 
due to Breit-Rabi-Zeeman splitting with 


7 = 1 


iHF 


f gsBo \ 

\hdvhfj 


( 10 ) 
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To find the change in 62 in the presence of the quadratic 
Zeeman effect one simply replaces the original (Zeeman- 
free) scattering length a in ([ 1 ]) with its renormalized form 
a* given by m, which results in a nonlinear contribution 
to the shift of the critical temperature, specifically 


at, 

rp 


~ 62 



( 11 ) 


Furthermore, by using , one can rewrite (HU in terms 
of the original scattering length a and renormalized am¬ 
plitude 62 s-s follows 


at, 

TO 


~ h* 
— 



( 12 ) 


where the coefficient due to the Breit-Rabi-Zeeman con¬ 
tribution is 


b 


* 

2 


l^b2 


(13) 


with 7 defined earlier. 

Let us consider the particular case of the Bq ~ 403G 
resonance of the hyperfine \F,mF >= | 1,1 > state of 
For this case ^|, S = ^12, mp = I = 3/2, and 
dPhf — 4:68MHz. These parameters produce Ahf = 3/4 
and 7 ~ 1.5 which readily leads to the following estimate 
of the quadratic amplitude contribution due to the HF 
mediated Breit-Rabi-Zeeman effect, 62 — 2.2562 — 42.3 
(using the mean-field value 62 ~ 18.8 [3l|)) in a good 
agreement with the observations 32|. It is interesting to 
point out that the obtained value of 7 for BEG is a 
result of a practically perfect match between the two par¬ 
ticipating energies: Zeeman contribution at the Feshbach 
resonance field, fJ-sBo ~ 4 x 10“^® J, and the contribu¬ 
tion due to Breit-Rabi hyperfine splitting between two 
ground states, hSvhf — 3 x 10“^® J. 

And finally, an important comment is in order regard¬ 
ing the applicability of the present approach (based on 
the Taylor expansion of (2)) to the field-induced modi¬ 
fication of the linear contribution (defined via the am¬ 
plitude 61 in (1)) to the shift in T,. According to the 
experimental curve depicting AT, vs a/Xp behavior, the 
linear contribution is limited by 10“^ < o/At < 5 x 10“^. 
Within the Feshbach-resonance interpretation, this cor¬ 
responds to a low-field ratio a/a^g ~ 1 which invalidates 
the Taylor expansion scenario based on using a small pa¬ 
rameter abgja <C 1 applicable in high fields only. Besides, 
as we have demonstrated earlier, the linear Zeeman ef¬ 
fect (valid at low fields only) is not responsible for any 
tangible changes of BEG properties. Therefore, another 
approach is needed to properly address the field-induced 
variation (if any) of the linear contribution 61 . 

To conclude, it was shown that accounting for 
a hyperfine-interaction induced Breit-Rabi nonlinear 
(quadratic) Zeeman term in the mean-field approxima¬ 
tion can explain the experimentally observed shift in the 


critical temperature T, for the condensate. It would 
be interesting to subject the predicted universal relation 
(13) to a further experimental test to verify whether or 
not it can also explain the shift in other bosonic-atom 
condensates. 
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